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ABSTRACT 

We explore the minimal conditions which enable the formation of metal-enriched 
solar and sub-solar mass stars. Using a one-zone semi-analytical model, we accurately 
follow the chemical and thermal evolution of the gas with the aim of understanding 
how the initial metal and dust content alters the cooling and fragmentation properties, 
hence the characteristic stellar mass. We find that in the absence of dust grains, gas 
fragmentation occurs at densities nn ~ [10 4 — 10 5 ] cm -3 when the metallicity exceeds 
Z ~ 10~ 4 Z Q . The resulting fragmentation masses are ^ 10M©. The inclusion of Fe 
and Si cooling does not affect the thermal evolution as this is dominated by molecular 
(mostly OH, H2O and CO) cooling even for metallicities as large as Z = 10~ 2 Zq. 
The presence of dust is the key driver for the formation of low mass stars. We focus 
on three representative core-collapse supernova (SN) progenitors (a Z — star with 
2OM and two Z = 10~ 4 Z Q stars with 20M© and 35M©), and consider the effects 
of reverse shocks of increasing strength: these reduce the depletion factors, /dc P = 
-Mdust/ (-^dust + M met ), alter the shape of the grain size distribution function and 
modify the relative abundances of grain species and of metal species in the gas phase. 
We find that the lowest metallicity at which fragmentation occurs is Z = 10~ 6 Zq for 
gas pre-enriched by the explosion of a 20M Q primordial SN (/d ep ^ 0.22) and/or by 
a 35M Q , Z = 1Q~ 4 Zq SN (/d ep ^ 0.26); it is ~ 1 dex larger, when the gas is pre- 
enriched by a Z — 10~ 4 Zq, 20 M Q SN (/d ep ^ 0.04). Cloud fragmentation depends on 
the depletion factor and it is suppressed when the reverse shock leads to a too large 
destruction of dust grains. These features are all consistent with the existence of a 
minimum dust-to-gas ratio, P cr , above which fragmentation is activated. We derive a 
simple analytic expression for P cr which depends on the total grain cross-section per 
unit mass of dust; for grain composition and properties explored in the present study, 
T> cr = [2.6 — 6.3] x 10~ 9 . When the dust-to-gas ratio of star forming clouds exceeds 
this value, the fragmentation masses range between O.OIMq and 1M , thus enabling 
the formation of the first low-mass stars. 



Key words: Cosmology: galaxies: evolution, stars: formation, Population II, ISM: 
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1 INTRODUCTION 

The formation of the first low-mass and long-lived stars 
marks an important step in cosmic evolution. Despite that 
a large fraction of the Galactic Halo has now been surveyed 
(Beers & Christlieb 2005; Helmi 2008), no evidence for the 
existence of primordial low-mass stars has been found to 
date. The most primitive star of the Galactic Halo has a 
[Fe/H] = -4.99 and, contrary to other extremely iron-poor 
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stars, shows no carbon, nitrogen or oxygen enhancement, 
resulting in a total metallicity of Z = 4.5 x 10~ Zq (Caffau 
et al. 2011). 

This lends support to the idea that the formation of 
solar or sub-solar mass starfl requires some minimal condi- 
tions which do not rely on carbon and oxygen line-cooling 

For WMAP7 cosmological parameters, the current Hubble time 
is tft = 13.75 Gyr (Jarosik et al. 2011). Depending on the stellar 
initial metallicity, this corresponds to the lifetime of stars with 
masses in the range [0.7 — 0.9]Mq (Raiteri et al. 1996) 
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(Frebel, Johnson, and Bromm 2007) and that can operate 
at very low but non-zero metallicity (Schneider et al. 2003, 
2006). 

This empirical evidence is supported by a solid theoreti- 
cal argument: the micro-physics of the H2 molecule, the only 
available coolant in ~ 10 6 Mq mini-halos at 2 ~ 20 — 30, sets 
a characteristic thermodynamic state where gravitational 
fragmentation of star-forming clouds is halted, which corre- 
sponds to a minimum Jeans mass of ~ 500Mq (see Bromm 
et al. 2009 and references therein) . Simulations starting from 
cosmological initial conditions have confirmed that vigorous 
fragmentation of collapsing clouds does not occur (Abel, 
Bryan & Norman 2002; Bromm, Coppi & Larson 2002; 
Yoshida et al. 2006; Yoshida, Omukai & Hernquist 2008), 
although there have been evidences for the formation of mas- 
sive binary systems (Turk, Abel, O'Shea 2009). 

The subsequent phase of gas accretion onto the cen- 
tral proto-stellar core is subject to radiative feedback ef- 
fects from the growing protostar which leads to a plausible 
but uncertain range of final stellar masses, [30 — 100] Mq 
(Omukai & Palla 2003; Tan & McKee 2004; Bromm & Loeb 
2004; Yoshida et al. 2008; McKee & Tan 2008; Hosokawa & 
Omukai 2009; Ohkubo et al. 2009). The accretion disk may 
become susceptible to fragmentation, forming groups of pro- 
tostars (Machida et al. 2008; Stacy, Greif & Bromm 2010; 
Clark et al. 2011; Greif et al. 2011), but current simulations 
are not able to follow the subsequent interactions of group 
members with the surrounding gas cloud. 

Rotation, turbulence and magnetic fields are known to 
play an important role in present-day star formation and 
their importance in the first star forming regions at high 
redshift still needs to be fully assessed (see e. g. Schleicher et 
al. 2009; Clark et al. 2011b and references therein). Whether 
or not these additional physical processes might enable the 
formation of low-mass Pop III stars remains an open ques- 
tion. 

Here we follow a different approach: guided by the 
observed constancy of the stellar characteristic mass in 
the local Universe, M c h ~ 0.3M Q (for recent reviews see 
Elmegreen 2009; Bastian, Covey & Meyer 2010), which has 
been interpreted to originate from the minimum Jeans mass 
imprinted during gravitational fragmentation (Li, Klessen & 
McLow 2003; Larson 2005; Jappsen et al. 2005; Elmegreen, 
Klessen & Wilson 2008), we explore the most extreme en- 
vironmental conditions where low-mass star formation can 
occur in the first mini-halos at high redshifts. 

The transition in the characteristic stellar mass, often 
referred to as Pop III/II transition, has been explored by 
many authors. On the relatively small scale of a single star 
forming cloud, semi-analytic and numerical studies have em- 
phasized the role of metals and dust. 

The effects of metal-line cooling, predominantly by 
O, C, Si and Fe, are relevant at relatively low densities, 
< [10 4 - 10 5 ] cm" 3 , before and around the onset of NLTE- 
LTE transition of the corresponding level populations. When 
the initial metallicity of the star forming cloud exceeds a 
critical value of Z CI ~ [10~ 4 - 1O" 3 ]Z , the minimum Jeans 
mass achieved by gravitational fragmentation is found to 
vary in the range [10 — 100] Mq, depending on the pres- 
ence/absence of molecular coolants in the chemical network 
and on the number of heating/cooling processes included in 
the modelling (Omukai 2000; Bromm et al. 2001; Schnei- 



der et al. 2002; Bromm & Loeb 2003; Santoro & Shull 
2006; Jappsen et al. 2007, 2009; Smith, Sigurdsson & Abel 
2008; Hocuk & Spaans 2010). In this case, the formation 
of sub-solar or solar mass stars, which may survive to the 
present time, is determined during the subsequent evolu- 
tionary phases, by gas accretion within each clump and dy- 
namical interactions within the cluster (Machida et al. 2009; 
Omukai, Hosokawa & Yoshida 2010). 

As in present-day conditions, thermal emission by colli- 
sionally excited dust grains is an efficient cooling channel at 
high densities, > [10 12 - 10 14 ] cm . As a result, the typical 
Jeans masses which can be achieved by gravitational frag- 
mentation are much smaller, [0.01 — 1] Mq (Schneider et 
al. 2003, 2006; Omukai et al. 2005; Tsuribe & Omukai 2006; 
Clark, Glover & Klessen 2008; Omukai et al. 2010; Dopcke et 
al. 2011). Thus, efficient dust-cooling leads to a transition in 
the characteristic mass by several orders of magnitude with 
respect to the primordial case. The metallicity at which this 
transition takes place depends on the dust depletion factor, 
defined as, 

/dcp ~ M dust + M mct ' W 

where M dust and M mot are the dust and gas-phase metals 
masses, respectively. 

In the model presented in Omukai et al. (2005), where 
we assumed present-day ISM composition and grain prop- 
erties, /dep = 0.49 (Pollack et al. 1994) and the onset of 
dust-induced cooling and fragmentation leading to the for- 
mation of ~ lM a fragments starts to be effective when 
10~ 6 Z Q < Z CI < 1O~ 5 Z . This result has been confirmed by 
numerical simulations (Tsuribe & Omukai 2006, 2008; Clark 
et al. 2008; Omukai et al. 2010; Dopcke et al. 2011). 

Since the properties of dust grains at high redshifts 
might be different with respect to those observed locally, 
in Schneider et al. (2006) we have investigated the dust- 
driven transition under the hypothesis that the collapsing 
star-forming cloud had been pre-enriched by metal-free core- 
collapse supernova (with a progenitor of 22Mq) and by a 
pair- instability supernova (with a progenitor of 195M©). 
The adopted metal composition and dust-grain properties 
were computed according to the model of dust nucleation 
in supernova ejecta developed by Todini & Ferrara (2001) 
and by Schneider, Salvaterra & Ferrara (2003). The deple- 
tion factors were /d ep = 0.23 (for the 22M Q progenitor) and 
/dc P = 0.65 (for the 195Mq) and low-mass fragmentation 
was found to occur already at Z cr = 10~ 6 Zq. 

An important confirmation of a Pop III/II transition 
driven by dust cooling has been recently provided by Caffau 
et al. (2011), who have reported the observation of a Galactic 
Halo star with Z = 4.5 x 10~ 5 Zq. Since this paper has 
appeared at the time of submission of the present work, we 
defer a specific analysis of this particular star to a future 
study. 

Here we intend to investigate the effects of varying the 
initial /d op on the gravitational collapse and fragmentation 
properties of pre-enriched star forming clouds. The ques- 
tions we intend to address are: (i) what is the minimum 
fragment mass that can be achieved by metal-line cooling? 
and (ii) What is the minimum dust depletion factor to ac- 
tivate dust-induced cooling? To answer these questions, we 
first calculate the thermal evolution of collapsing gas clouds 
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that cool only through molecular (H2, HD, H2O, CO, OH) 
and metal (CI, CII, 01, Sill, Fell) line cooling. Then, we 
use the model developed by Bianchi & Schneider (2007) to 
compute the mass, composition and size distribution of dust 
formed in the ejecta of a few representative SN progenitors. 
Starting from these initial conditions, we take into account 
the partial destruction of the newly formed dust by the SN 
reverse shock. In particular, we consider reverse shocks of 
increasing strength, assuming increasingly large circumstel- 
lar densities. Following the models by Bianchi & Schneider 
(2007) we can control the reduction of the depletion factor 
(more than 80% of the initial dust mass is destroyed) in a 
self-consistent way and explore the implications for the frag- 
mentation properties of pre-enriched star forming clouds. 

The plan of the paper is as follows: in Section 2 we 
present the model; in Section 3 we discuss the thermal evo- 
lution driven by line-cooling and the resulting fragmenta- 
tion properties. In Section 4 we describe the adopted model 
grid of core-collapse SN progenitors and the associated dust 
properties. In Section 5 we present the thermal evolution and 
fragmentation of clouds pre-enriched by these explosions. In 
Section 6 we discuss the results and, finally, in section 7 we 
summarize our main conclusions. 



2 THE MODEL 

The thermal and chemical evolution of collapsing gas clouds 
is investigated using the model presented in Omukai (2000) 
and in Omukai et al. (2005) with CO, H2O and OH cooling 
rates recently modified as in Omukai et al. (2010). As a new 
feature of the code, we have included the contribution of Fe 
and Si to metal-line cooling. Here we briefly summarize the 
features of the model which are relevant for the purposes of 
the present study and we refer the interested reader to the 
original papers for a thorough description of the code. 

The temperature evolution is computed by solving the 
energy equation, 



de 
dt 



d 1 
dt p 



A ne 



(2) 



where the pressure, p, and the specific thermal energy, e, are 
given by, 



pkT 



where / p is the ratio between pressure gradient to gravity 
(see Omukai et al. 2005 for the expression) and the net cool- 
ing rate, A not , which consists of three components, 



A ne t — Ali nc + A CO nt + A c 



(6) 



The first component, Aii nc , represents the cooling rate due 
to the emission of line radiation, which includes molecular 
line emission of H2, HD, OH, H2O and CO (Omukai et al. 
2010) and atomic fine-structure line emission. In addition 
to CI, CII and OI (Omukai et al. 2005), here we have con- 
sidered the contribution of Fe and Si lines (see section [3] 
below). The second component, A cont , represents the cool- 
ing rate due to the emission of continuum radiation, which 
includes thermal emission by dust grains and H2 collision- 
induced emission. The last term, A cnem , indicates the cool- 
ing/heating rate due to chemical reactions. When the gas 
cloud is optically thick to a specific cooling radiation, the 
cooling rate is correspondingly reduced by multiplying by 
the photon escape probability (see Omukai 2000). 

We do not consider the effects of the Cosmic Microwave 
Background (CMB) radiation field, which provides a redshift 
dependent effective temperature floor (Clarke & Bromm 
2003; Omukai et al. 2005). At z ^ 10 - 15, the evolution 
of star-forming clouds can be significantly affected but only 
at relatively large metallicities, Z ^ 10~ 2 Z@ (Smith et al. 
2009; Schneider & Omukai 2010). We will return to this 
point in section 6. 

The effects of dust grains on the thermal evolution are 
due to (i) dust thermal emission and (11) the increased for- 
mation rate of H2 on its surface and (Hi) the heating asso- 
ciated to the latter process. The treatment of these physical 
processes is the same as in Schneider et al. (2006). Here we 
report only the fundamental equations which regulate dust 
thermal emission and refer the interested reader to Schnei- 
der et al. (2006) for further details. 



2.1 Dust cooling 

Dust grains affect the thermal evolution of the collapsing gas 
cloud when the energy transfer between gas and dust due to 
collisions is radiated away. Dust thermal emission rate can 
be expressed as, 



(3) 



A. 



grL erg/s/gr] 



4aT k p (3 C: 



Pgr 



(7) 



kT 



7 — 1 p m H 



(4) 



and p is the central density, T is the temperature, 7 is the 
adiabatic exponent, p is the mean molecular weight and 
ma is the mass of hydrogen nuclei. 

The terms on the right-hand side of the energy equation 
are the compressional heating rate, 



dp = P 
dt t co \ 



with t co i = 



VWpV 32Gp 



3tt 



(5) 



2 The mean molecular weight is 1.23 for a fully atomic gas and 
2.29 for a fully molecular gas. 



where a is the Stefan-Boltzmann constant, T gr is the dust 
temperature, Kp is the Planck mean opacity of dust grains 
per unit dust mass, /3 C sc is the photon escape probability, 
and the dust mass density, p gT , can be expressed as pT>, 
with T> indicating the dust-to-gas ratio. Following Omukai 
(2000), the photon escape probability can be written as, 



(3 CBC = min(l, r ) 



with 



r = KppXj, 



(8) 



where the optical depth is estimated across one local Jeans 
length (roughly the size of the central core region in the 
Penston-Larson similarity solution) . The Planck mean opac- 
ity per unit dust mass is, 



Kp = 



B v (T ST )K v dv 



(9) 
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Figure 1. Thermal evolution of collapsing clouds with initial 
metallicities Z = 0, 10~ 4 Z Q , 10 _3 Z Q , 1O _2 Z , from top to bot- 
tom. For each of the Z > models, two lines are shown: the 
thin one represents the prediction of the Omukai et al. (2005) 
model; the thick one shows the effects of including [Fell] and 
[Sill] cooling. The diagonal dotted lines indicate locii of constant 
Jeans masses, as labelled in the figure, and the points with trian- 
gles mark the states where fagmentation conditions are met (see 
text). 



where k„ is the absorption coefficient for the frequency v 
per unit dust mass and B v (T gI ) is the Planck brightness 
(see Schneider et al. 2006). 

The dust temperature, T gI , is determined by the balance 
between dust thermal emission and the dust heating rate due 
to collisions with gas particles (Hollenbach & McKee 1979), 
H gr — pA gr , with 

H gr [erg/s/cm 3 ] = ^kT~2kT gI ) ^ (1Q) 

Icoll 

where t~ oll = riH a gr vh f is the average time between two 
successive collisions, n gr and a gI are the grain number den- 
sity and cross section, vn is the average speed of hydrogen 
nuclei and / = nv/nnVH ~ 0.395 measures the contribution 
of other specie^. 

In the density range where collisional dust heating is rel- 
evant (riH > 10 12_13 cm -3 ), it is a good approximation to as- 
sume that the gas is fully molecular and follows a Maxwellian 
distribution so that, 




If we further write, 

n S v o-g r = n H m H (l + 4j/ He ) ST> (12) 

in eq. (|10[) . where S is the total grain cross section per unit 
mass of dust, the equation pA gr = H gI appears to be inde- 
pendent of the dust-to-gas ratio. Thus, T gr depends only on 
gas temperature and density and on grain specific properties 
(cross section and mean opacity, see Schneider et al. 2006 
for further details). 




Log (n H /cm- 3 ) 



Figure 2. Individual metal fine-structure lines (Fell, Sill, 
CII+CI, OI) and molecular lines (H 2 , CO, OH+H 2 0) contribu- 
tions to the cooling rate of a gas cloud with initial metallicity 
Z = 10~ 2 Zq (see FigQJ. 



3 FRAGMENTATION BY METAL-COOLING: 
EFFECTS OF FE AND SI LINES 

In this section, we discuss the thermal evolution of collapsing 
gas clouds when only molecules and metals contribute to 
the cooling rate in eq. ((S}. Thus, we assume that no dust 
is present. For simplicity, we adopt an ISM with the solar 
elemental composition (Anders & Grevesse 1989) and the 
mass of each element is decreased proportionally to the total 
metallicity. 

To estimate the effects of Fell and Sill lines on gas cool- 
ing, we do not solve the full Fe and Si chemistry (as we do 
for the other elements) but we simply assume that all the Fe 
and Si contributes to cooling. Fig.[T]shows the thermal evolu- 
tion of collapsing gas clouds with varying initial metallicities 
(Z = 0, 10~ 4 Z q ,10- 3 Zq,W- 2 Z q , from top to bottom). For 
the Z > models, thin lines represent the predictions ob- 
tained with Omukai et al. (2005) models and thick lines show 
the current models with the inclusion of Sill and Fell line- 
cooling. As it is apparent from the figure, the difference is 
very small and confined to the density range nu < 10 6 cm -3 . 
In fact, even for the highest metallicity model that we have 
considered, Z — 10~ 2 Zq, the cooling rate is dominated by 
molecular lines, specifically by H2O and OH. This is shown 

3 Since we are interested in very low gas metallicities (10 _7 Zq ^ 
Z ^ 10~ 4 Zq) we can safely neglect the contribution of heavy 
elements and write / = 0.3536 + 0.5j/He with yn e = nne/ n n = 
0.083, which corresponds to a He mass fraction of Y=0.25. 
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Figure 3. Mass of dust predicted by SN as a function of the progenitor stellar mass for two initial metallicities: Z = (left panel) and 
Z = 1O~ 4 Z0 (right panel). In each panel, the solid line refers to the mass of dust predicted at the end of the condensation phase (norev) 
and for the revl, rev2, rev3 cases (dotted, dashed, long-dashed lines, respectively) correponding to ISM densities of 10~ 25 , 10 -24 and 
10 — 23 gr/cm — 3 , respectively. 



Table 1. Ejecta composition of the three SN models; for each of these we show the depletion factor, /d Gp (first column), the mass of 
dust grains, M;j ust (second column), the total mass of gas-phase metals, M mc t (third column) and the mass of O and C atoms (fourth 
and fifth columns, respectively) in the norev, revl, rev2, and rev3 cases. All masses are in solar units. The last column shows the 
total dust cross section per unit dust mass (see section [2. 111 . 
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in Fig. [2] where we compare the individual contributions to 
the cooling rate. Although Fell and Sill are the most im- 
portant coolants among metal fine-structure lines in some 
density ranges, their contribution to the total cooling rate 
is smaller than that provided by molecules. 

In Fig[T] the diagonal dotted lines represent constant 
Jeans mass values in the (riH, T) plane. These are computed 
using the relation, 



il/.T 



1.69 x 10 M© 



(1.23) 



- 3 / 2 / T 
V200K 



J \ 10 4 cm 3 J 



(13) 

-1/2 



where the Jeans length is, 

A = / 7rfc B r \ 1/2 

\Gfj,m K p J 



(14) 



The points with triangles mark the states where fragmen- 
tation conditions are met. Following Schneider & Omukai 
(2010), we identify the fragmentation epochs by requiring 
that the adiabatic index becomes 7 > 0.97 after a phase of 
cooling (where 7 < 1). To eliminate "false" fragmentation 
points, where 7 is less than 1 only for a short period of time, 
we impose that 7 < 0.8 during the cooling phase preceding 
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Figure 4. Size distribution function of the dominant dust species resulting from the explosion of a Z = 0, Af s tar = 20-Mq (left panel), 
Z = 10 — 4 Zq, Mstar = 2OM (central panel), and Z = 1O _4 Z , M a t ar = 35M (right panel). Solid lines show the distribution at the end 
of the condensation phase (norev) and dotted, dashed, long-dashed lines illustrate the effect of reverse shocks with increasing strength 
(revl, rev2, rev3, respectively). 



the fragmentation epoch (see Schneider & Omukai 2010 for 
more details). 

It is clear from the figure that the minimum fragmen- 
tation masses that can be achieved by metal and molecular 
line-cooling are ^ IOMq, even for gas clouds with initial 
metallicity as high as Z = W~ 2 Zq. Thus, the additional 
contribution of Fe and Si lines to the chemical network and 
cooling rate do not change the main conclusions of previous 
investigations, based on the model presented in Omukai et 
al. (2005). In particular, the formation of solar or sub-solar 
mass fragments requires the presence of dust grains and 
dust-induced cooling at densities tih > 10 12 cm -3 . 



4 MODEL GRID OF CORE-COLLAPSE SN 

Figure shows the masses of dust produced by core col- 
lapse SN as a function of the progenitor star mass and initial 
metallicities equal to Z = (left panel) and Z = 1O _4 Z 
(right panel). The solid lines show dust masses at the end of 
the condensation phase (no reverse shock model) according 
to the model developed by Todini & Ferrara (2001) and re- 
visited by Bianchi & Schneider (2007). Here we refer to the 
fiducial case, where we assume that the first seed critical 
dust clusters are formed with a minimum of two monomers 
and that the sticking coefficient is equal to 1 (we refer to 
Bianchi & Schneider 2007 for more details). We see that the 
most massive SN (M ata r > 25Mq) are less efficient sources 
of dust due to the strong fallback of metals during the explo- 
sion. This is particularly true for the primordial case, where 
metal (and therefore dust) production is zero beyond a pro- 
genitor stellar mass of 3OM . The dotted, dashed, and long- 
dashed lines show the masses of dust after the passage of the 
reverse shock when the surrounding ISM density is assumed 
to be equal to 10~ 25 , 10~ 24 and 10~ 23 gr/cm -3 , respectively. 
The reverse shock is modelled following the prescription of 
Bianchi & Schneider (2007) and dust destruction is caused 
by both thermal and non-thermal sputtering (see the origi- 
nal paper for more details). The fraction of the initial dust 
mass which is destroyed and returns back to the gas phase 



increases with the ISM density. Hereafter, we will refer to 
the ejecta properties at the end of the condensation phase 
as the norev case; with revl, rev2, rev3 cases we will in- 
dicate the models with reverse shock of increasing strength 
(increasing ISM density). 

As it is reasonable to expect, the mass of dust re- 
leased by SN increases with the metallicity of the progen- 
itor stars. We focus here on three representative SN mod- 
els (highlighted in the figure): a Z — 0, 20 Mq progenitor 
(2OM , 0) and two Z — 10 -4 Zq progenitors with masses 
of 20 (2OM , W- i Z e ) and 35M (35M , 1O _4 Z ). Ta- 
ble [1] describes the final ejecta composition for these three 
SN in the norev, revl, rev2 and rev3 cases. The first 
column lists the depletion factor, /de P , with values in the 
range 4.92 x 10~ 3 ^ /d cp 0.91, decreasing with reverse 
shock strength. Thus, the conditions we will explore span 
the value appropriated for the local ISM (/dcp = 0.49) as 
well as values which are 2 dex smaller. Note that the deple- 
tion factors are systematically larger for the (2OM , 0) and 
the (35M , KT 4 Z ) SN progenitors, despite the smaller fi- 
nal dust and metal yields with respect to the SN progenitor 
(2OM , 1O" 4 Z ). 

The elemental composition of the ejecta depends on the 
SN progenitor and reverse shock model. As it can be seen 
from Table [TJ oxygen is always the dominant element at the 
end of the condensation phase (norev cases) and carbon is 
returned to the gas phase after the partial destruction of 
graphite grains by the reverse shock (revl-rev3 cases). We 
therefore expect that the thermal properties of gas enriched 
by these ejecta will depend on /d ep due to variations both 
in the abundance of dust grains and in the composition of 
gas-phase metals. 

The reverse shock significantly affects the grain size dis- 
tribution function, decreasing the number of larger grains 
and flattening the distribution towards smaller grain sizes. 
This is shown in Fig(4] where, for each model, the domi- 
nant grain species have been reported. Only amorphous car- 
bon (AC) grains form in the ejecta of the (2OM ,O) and 
the (35M ,1O _4 Z ) progenitors; additional grain species 
form in the ejecta of the (2OM ,1O _4 Z ) progenitor, among 
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Figure 5. Planck mean opacity per unit dust mass as a func- 
tion of grain temperature. The solid line refers to models Z = 0, 
M star = 2OM , and Z = 1O~ 4 Z , M star = 35M which have 
very similar dust properties. Also, for these two models we only 
show the opacity at the end of the condensation phase because 
the passage of the reverse shock has a negligible effect. The 
dashed region illustrates the maximum (norev) and minimum 
(rev3) Planck mean opacity per unit mass for the Z = 1O _4 Z , 
M a tar = 2OM model. Line breaks are due to grain sublimation 
(see text). 

which FesCU, MgSi03,SiC>2, and AI2O3 (see also Bianchi & 
Schneider 2007 for more details). Thus, the former two SN 
progenitors share very similar dust properties. This is re- 
flected by their corresponding grain cross sections per unit 
dust mass (see the last column in Table [1} and Planck mean 
opacities (see eq. [5| which are plotted in Fig. [5] as a func- 
tion of dust grain temperature. The solid lines show the 
predicted opacity for (20A/ O ,0) and (35M , 10~ 4 Z Q ) pro- 
genitors with no variations induced by dust destruction in 
the reverse shock; the shaded region indicates the maximum 
(norev case) and minimum (rev3 case) Planck mean opac- 
ity predicted for the (2OM , 1O _4 z7 ) progenitor. The line 
breaks are due to grain sublimation, which occurs at a tem- 
perature of 1810 K for AC grains and of 1140 K for the ad- 
ditional grain species (Kozasa, Hasegawa, & Nomoto 1989). 
Since the grain cross section and Planck mean opacity con- 
trol the efficiency of dust cooling (see Section \2. II) . in what 
follows we restrict the discussion to the most dissimilar SN 
progenitors, hence a 2OM star with metallicities Z = and 
Z = 1(T 4 Z©. 



5 RESULTS 

In this section we investigate the thermal evolution and frag- 
mentation of collapsing gas clouds enriched by each of the 
two representative supernova progenitors. We assume that 
their ejecta uniformly enrich nearby collapsing gas clouds. 



Depending on its total gas mass, the collapsing cloud can 
be enriched to different total metallicities, where by total 
metallicity we mean the total mass of metals (in the gas- 
phase and in dust grains) relative to the mass of gas. Figs[6] 
and show the thermal evolution of a collapsing gas cloud 
with initial metallicity equal to 10 -7 Zq, W~ 6 Zq, 1O _5 Z , 
and W~ 4 Zq (from top left to bottom right); in each panel, 
we plot different reverse shock models varying the deple- 
tion factor and the ejecta composition self-consistently (from 
bottom to top norev, revl, rev2, rev3 correspond to solid, 
dotted, short-dashed and long-dashed lines); as a reference 
case, the evolution predicted for a Z = gas is also shown 
(top solid line). 

The comparison among different panels shows that the 
deviation from the primordial gas evolution increases with 
metallicity; when Z — 1O~ 7 Z , this is limited to rela- 
tively large densities, nu ^ 10 12 cm -3 , but progressively 
extends to lower densities as the metallicity increases. At 
metallicities 1Q~ 7 Zq ^ Z ^ 1O _4 Z the evolution is dom- 
inated by dust through both direct cooling (high densities; 
> 10 11 cm~ 3 ) and increased H2 cooling due to H2 formation 
on dust grains (lower densities; < 10 9 cm -3 ), as already dis- 
cussed in Omukai et al. (2005) and Schneider et al. (2006). 
At each given metallicity the largest deviation from the pri- 
mordial case is associated to the norev case, characterized 
by the largest /d cp , and the evolution progressively tend to 
the primordial limit when lower depletion factors are con- 
sidered (models revl, rev2 and rev3). 

The points with triangles in Figs[6]and[7]mark the states 
where fragmentation conditions are met. The lowest metal- 
licity at which fragmentation occurs is Z — 1O _6 Z for the 
(2OM , 0) SN progenitor and it is 1 dex larger, Z = 1O~ 5 Z , 
for the (2OM0, 1O _4 Z ) SN progenitor. In both cases, how- 
ever, cloud fragmentation depends on the depletion factor 
and it is suppressed when the strength of the reverse shock 
leads to a too large destruction of dust grains (rev2 and 
rev3 cases). 

These features are all consistent with the existence of 
a minimum dust-to-gas ratio above which fragmentation is 
activated. Given our definitions of metallicity, Z, and dust 
depletion factor, /d ep , the dust-to-gas ratio can be written as 
T> — Zfdcp- Thus, the minimum T> is reached at a lower total 
metallicity for models with larger depletion factors, such as 
the (20M Q ,0) SN progenitor (see Table [T]). 

FiglU represents the results of all the models investi- 
gated so far in a Z - /d ep plane. Empty symbols indi- 
cate no fragmentation and filled symbols represent mod- 
els where fragmentation conditions are met. The solid line 
marks the border between the no-fragmentation and frag- 
mentation zones and corresponds to a minimum dust-to-gas 
ratio of only T> CT — 4.4 x 10 -9 . At any given metallicity, 
low-mass fragmentation can be activated when T> D cr or 

/dcp ^ "Z^cr jZ . 

A simple argument for the origin of T> cr runs as fol- 
lows. The second high-density dip in the thermal evolution 
of collapsing gas clouds is activated when the energy transfer 
rate between gas and dust is more efficient than the com- 
pressional heating rate due to gravitational contraction, that 
is, 

H v /p>-p±\. (15) 
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Figure 6. Thermal evolution of clouds enriched by the (2OAf0,O) SN model. In each panel, we have assumed that the gas has been 
pre-enriched to a total metallicity of 10 _7 Zq, 10 — ®Zq, 10~ 5 Zq, and 10 — 4 Zq (from top left to bottom right). In each panel, we show 
the evolution for four values of /d ep corresponding to the norev, revl, rev2, and rev3 cases (bottom solid, dotted, short-dashed, 
long-dashed, respectively). For each of these cases, the ejecta composition is varied consistently, as presented in Table [T] and Figs. [4] ED 
As a reference model, we also show the evolution predicted for a purely primordial gas cloud (top solid line). The diagonal dashed lines 
represent lines of constant Jeans mass and the triangles identify the points where fragmentation conditions are met (see text). 



Using eqs. (fT0|) - (fT2)) for the left-hand term and eqs. ©-([SI) 
for the compressional heating rate, eq.(|15[) leads to a condi- 
tion on the dust-to-gas ratio times the cross section per unit 
dust mass, 



SV > 1.4 x 10~ 3 cm 2 /gr 



- T - 


-1/2 


nH 


-1/2 


Ll0 3 K_ 




_10 12 cm- 3 . 





-(16) 



where we have assumed that at the onset of dust cooling, 
before gas and dust become thermally coupled, T gr << T. 
The above relation has been evaluated at reference gas tem- 



perature and density of T ~ 10 3 K and nu ~ 10 12 cm -3 . 
Indeed, the sharp temperature drop, which activates frag- 
mentation, can be achieved only if dust cooling exceeds the 
compressional heating rate when tih < 10 12 cm -3 and H2 
formation heating keeps the gas temperature at ~ 10 3 K. 
Beyond this density, H2 becomes fully molecular, chemical 
heating is no longer effective, and dust-cooling causes only 
a mild temperature decrease. It follows from eq. (|16|l that 
the critical dust-to-gas ratio is inversely proportional to the 
total grain cross section per unit mass of dust. As it can 
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Figure 7. Same as Fig.[g]but for the (20Af Q , 10 _4 Z Q ) SN model. 



be inferred from the last column of Table [TJ differences in 
S among different reverse shock models for a given SN pro- 
genitor are less than 15%. Taking S to vary in the range 
2.22 < 5/10 5 cm 2 /gr < 5.37, the condition on the dust-to- 
gas ratio reads, 



9 


- j, - 


-1/2 


n H 


-1/2 




_10 3 K_ 




_10 12 cm" 3 _ 





consistent with value of T> cr inferred above and drawn in the 
figure. 

It is important to stress that this simple analytic ar- 
gument is not able to capture all the important details of 
the evolution which may activate dust-induced fragmenta- 
tion. In particular, the condition that 7 < 0.8 during the 



dust-induced cooling phase depends also on the grain opac- 
ity: if we increase the dust opacity for fixed surface area, the 
dust-gas coupling temperature becomes lower, leading to a 
deeper dip in the equation of state. 

In addition, the fragmentation criteria that we follow 
(7 > 0.97 after a phase of cooling where 7 < 0.8) have 
been established through a comparison with numerical stud- 
ies (Tsuribe & Omukai 2006; Omukai et al. 2010). It is clear 
that the process of fragmentation does not depend only on 
the gas thermal evolution, being influenced also by the de- 
tailed hydrodynamical properties of the clouds, such as the 
initial amplitude of fluctuations, rotation and turbulence. 
Thus, our study can provide a guideline for more sophisti- 
cated numerical simulations. 
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Figure 8. Fragmentation (filled symbols) and no-fragmentation 
(empty symbols) conditions for each considered SN progen- 
itors and reverse shock model. Red circles, blue squares, 
and green triangles refer to (2OM ,O), (20A/ Q , 1O~ 4 Z ), and 
(35Mq, 10 — 4 Zq) models, respectively. In addition, we also show 
the SN and PISN models analyzed in Schneider et al. (2006, cyan 
and magenta exagons) and the model which adopts local ISM 
conditions as in Omukai et al. (2005, black exagon). The black 
solid line shows the predicted behaviour for a constant dust-to-gas 
ratio of V CI = 4.4 X IO" 9 (see text). 

6 DISCUSSION 

Although we have considered only star-forming regions in 
the first mini-halos, metals and dust rapidly enrich more 
massive halos, with masses M ~ lO 8 Af0 at z ~ 10, the so- 
called first galaxies (Bromm & Yoshida 2011). When the gas 
is accreted onto these halos, it is shock-heated at the virial 
radius of the parent dark matter halo to virial temperatures 
T v i r ~ 10 4 K; then it cools via the atomic-hydrogen Lyman- 
a line, falling inward and attaining increasingly larger den- 
sities. In Schneider et al. (2002), we investigate the depen- 
dence of thermal evolution of collapsing gas clouds on the 
initial virial temperature (halo mass). For metal-free gas, 
the evolution with temperature of the ionization fraction 
and of the fraction of molecular hydrogen is independent 
of the halo virial temperature. If the gas is fully ionized, 
molecules form through nonequilibrium recombination, lead- 
ing to overall fractions that are much higher than in the ex- 
panding homogeneous universe (Uehara & Inutsuka 2000). 
As a consequence of enhanced H2 and HD fractions, the gas 
rapidly cools to the CMB temperature and ~ IOMq frag- 
ments are formed (Johnson & Bromm 2006). 

In the presence of a strong LW background, the gas 
in the first galaxies remains free of H2 (Omukai 2001), it 
avoids fragmentation and leads to the formation of a single 
super-massive star with mass as high as ~ [10 5 — 1O 6 ]M0, 
which then evolves into a super-massive black hole (Bromm 
& Loeb 2003b). Omukai, Schneider & Haiman (2008) find 
that, even in these systems, fragmentation is inevitable when 



the gas is pre-enriched above a critical metallicity, whose 
value is between 3 x 10~ 2 Z@ (in the absence of dust) and as 
low as 5 x 10~ 6 Zq (with a dust-to-gas mass ratio of about 
0.01Z/Z®). 

Recently, Dekel & Birnboim (2006) and Dekel et al. 
(2009) have proposed an alternative accretion mode, the 
so-called 'cold accretion': gas falls in the halos along cos- 
mic filaments and it is deposited much closer to the center of 
the galaxy, at much higher densities. In these conditions, the 
fragmentation properties of the gas are likely to be modified. 
Safranek-Shrader, Bromm & Milosavljevic (2010) find that 
H2 and HD dominate gas cooling at T < 10 3 K and metal- 
licity has no effect on fragmentation. In the presence of a 
strong LW background which photo-dissociates H2, metal- 
line cooling becomes important and, when the metallicity is 
Z ~ W~ i Z e , ~ 3M fragments are formed. While this ap- 
pears an interesting possibility to form low-mass stars, the 
final fragmentation mass depends sensitively on the adopted 
initial gas density. Thus, a larger set of initial conditions 
needs to be investigated to draw robust conclusions. 

Finally, it is important to mention the effects of the 
CMB radiation field on the thermal evolution of collaps- 
ing gas clouds, which have not been included in the present 
study. In general, the presence of the CMB radiation field 
is important at high redshifts and metallicities (Smith et al. 

2009) , inhibiting fragmentation. Schneider & Omukai (2010) 
find that in the absence of dust grains, the fragmentation 
mass scales at z > 10 are always 100s of Mq, indepen- 
dent of the gas initial metallicity. When Z < W~ 2 Zq, dust- 
cooling remains relatively insensitive to the presence of the 
CMB, and sub-solar mass fragments are formed at any red- 
shift. It is only when Z > W~ 2 Zq and z > 10 that the 
heating of dust grains by the CMB starts to inhibit the for- 
mation of sub-solar mass fragments (Schneider & Omukai 

2010) . Thus, the conclusions of the present study remain 
valid in all but these extreme conditions. 



7 CONCLUSIONS 

In this study, we have explored the conditions which enable 
the formation of solar and sub-solar mass fragments during 
the collapse of pre-enriched gas clouds within the first mini- 
halos at z ~ 20 — 30. 

By carefully analysing the thermal and chemical evo- 
lution of star-forming clouds, we find that in the absence 
of dust grains, gas fragmentation occurs at densities rin ~ 
[10 4 - 10 5 ] cm" 3 when the metallicity exceeds Z ~ W~ 4 Z Q . 
The resulting fragmentation masses, identified by the mini- 
mum Jeans mass at the fragmentation epoch, are IOMq. 
The inclusion of Fe and Si cooling (which can dominate 
metal-line emission in some density ranges) does not af- 
fect the thermal evolution as this is dominated by molecular 
(mostly OH, H2O and CO) cooling even for metallicities as 
large as Z = 10~ 2 Z Q . 

The presence of dust is the key driver for the forma- 
tion of solar-mass stars. We have explored the fragmentation 
properties of collapsing clouds pre-enriched with dust and 
metals synthesized by core-collapse SNe. Focusing on three 
representative SN progenitors (a Z — star with 20A1q 
and two Z — 10~ 4 Zq stars with 20M© and 35Mq), we 
compute the mass of dust and metals (hence the depletion 
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factor), grain composition and size distribution using the 
model by Bianchi & Schneider (2007). To determine the min- 
imum dust depletion factor required to activate dust-driven 
fragmentation, we start from the above initial conditions 
and consider the effects of SN reverse shocks of increasing 
strength: these reduce the depletion factors, alter the shape 
of the grain size distribution function and modify the rel- 
ative abundances of grain species and of metal species in 
the gas phase. Following this approach, we can control all 
these effects in a self-consistent way and explore the impli- 
cations for the fragmentation properties of pre-enriched star 
forming clouds. 

We find that the lowest metallicity at which fragmen- 
tation occurs is Z = 10 _6 Zq for gas pre-enriched by the 
explosion of a 20Mq primordial SN and/or by a 35Mq, 
Z — 10~ 4 Zq SN; it is ~ 1 dex larger, when the gas is pre- 
enriched by a Z = 10~ 4 Z Q , 20 M SN. In both cases, how- 
ever, cloud fragmentation depends on the depletion factor 
and it is suppressed when the reverse shock leads to a too 
large destruction of dust grains. 

These features are all consistent with the idea that there 
is a minimum dust-to-gas ratio, D cr , above which frag- 
mentation is activated. By requiring that the energy trans- 
fer rate between gas and dust is larger than the compres- 
sional heating rate, we derive a simple analytic expression 
for D cr which depends on the total grain cross-section per 
unit mass of dust; for grain composition and properties ex- 
plored in the present study, this critical dust-to-gas ratio 
is £> cr = [2.6 - 6.3] x 10~ 9 . When the dust-to-gas ratio of 
star forming clouds exceeds this value, the fragmentation 
masses range between O.OIM© and IMq, thus enabling the 
formation of the first low-mass stars. 

According to this criterium, the recently observed 
Galactic Halo star with a total metallicity of Z — 4.5 x 
W' 5 Z Q (Caffau et al. 2011) requires a minimum dust de- 
pletion of /dc P ^ (4 — 9) x 10~ 3 in the parent gas cloud, con- 
sistent with the values predicted for ordinary core-collapse 
SN investigated in the present study. 
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